n=7; N=32
mi=c(5,5,4,6,5,5,5)
Mi=c(15,15,14,17,15,15,15)

# also need vector of means (ybari) and variances (s2i)
# look at example for calculation means and variances
# example is "2stgcluster"